The Local and the Occupation Time of a Particle Diffusing in a Random Medium 
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We consider a particle moving in a one dimensional potential which has a symmetric determin- 
istic part and a quenched random part. We study analytically the probability distributions of the 
local time (spent by the particle around its mean value) and the occupation time (spent above its 
mean value) within an observation time window of size t. In the absence of quenched randomness, 
these distributions have three typical asymptotic behaviors depending on whether the deterministic 
potential is unstable, stable or flat. These asymptotic behaviors are shown to get drastically modi- 
fied when the random part of the potential is switched on leading to the loss of self- averaging and 
wide sample to sample fluctuations. 



PACS numbers: 05.40.-a, 46.65.+g, 02.50.-r, 05.70.Ln 

How many days, out of a total number of t days, does a 
tourist spend in a given place? If the tourist is a Brown- 
ian particle, this local time (per unit volume) n t (f) spent 
by it in an infinitesimal neighbourhood of a point f in 
space has been of interest to physicists and mathemati- 
cians for decades. If R(t') denotes the position of the 
random walker at time t', the local time is defined as 



n t (r) 



R(t') 



dt'. 



(1) 



Clearly J n t (r)df = t. The local time is a very use- 
ful quantity with a variety of important applications in 
fields ranging from physics to biology. For example, in 
the context of polymers in a solution, the local time nt(r) 
is proportional to the concentration of monomers at r in a 
polymer of length t and can be measured via light scatter- 
ing experiments. In "True" self-avoiding walk (TSAW), a 
simple model of a polymer that takes into account the ex- 
cluded volume effect, the local time already spent by the 
walker at a site decides the future rate of revisits to that 
site 0] . The distribution of local time (LTD) also plays 
an important role in the study of diffusion in porous rocks 
P,|) . Another important application of the local time is 
in bacterial chemotaxis where the phenomenon of 'tum- 
bling' is quantified by the local time spent by a bacteria 
at a point M. 

A related quantity, also of wide interest to both physi- 
cists jq] and mathematicians p], is the occupation time 
T t (D) = J D n t (r)dr spent by the walker in a given region 
D of space. Recently the study of the occupation time 
has seen a revival due to its newfound applications in 
the context of persistence in noncquilibrium statistical 
physics [m]. The dynamics in such systems is typically 
modelled by a stochastic process x(t) whose statistical 
properties provide important informations about the his- 
tory of evolution in these systems jj. A quantity that 
acts as a useful probe to this history dependence is the 
occupation time T t = J 8 [x(t')]dt' , the time spent by 
the process on the positive side within a window of size 



t ||]. For instance, if x(t) represents a one dimensional 
Brownian motion, the probability distribution of the oc- 
cupation time (OTD) is given by the celebrated 'arcsinc' 
law of Levy §, Prob[T t = T] = l/ir^T(t-T) which 
diverges at the end points T = and T — t indicating 
the 'stiffness' of the Brownian motion, i.e., a Brownian 
path starting at a positive (negative) point tends to re- 
main positive (negative). In the mathematics literature 
this result has been extended in several directions using 
modern tools of probability theory such as martingale 
calculus and excursion theory |l0) . The OTD has been 
studied in a variety of physical systems. For example, it 
has been proposed as a useful tool to analyse the mor- 
phological dynamics of interfaces jll[] . Recently the OTD 
has also been used to analyse the experimental data on 
the "on-off" fluorescence intermittency emitting from col- 
loidal (CdSe) semiconductor quantum dots |{l2| . Exotic 
properties such as a phase transition in the ergodicity of 
a process by tuning a parameter have also been exhibited 
by the OTD in the diffusion equation 

While the LTD and the OTD have been studied ex- 
tensively for pure systems, they have received much less 
attention in systems with quenched disorder, the subject 
of this Letter. The study of these quantities can provide 
valuable informations in disordered systems. For exam- 
ple, if one launches a tracer particle in a system with 
localized impurities, the particle diffuses in this random 
medium, occasionally gets pinned in the region near the 
impurities till the thermal fluctuations lift it out of the 
local potential well and then it diffuses again. Since the 
local time spent by the particle at a given point in space 
is related to the concentration of impurities there, it can 
be used as a valuable probe to image the inhomogeneitics 
in a given sample. In this Letter we present exact results 
for the LTD and the OTD in the Sinai type models @ 
where a particle diffuses in a one dimensional random 
potential. Both of these distributions turn out to have 
very rich behaviors despite the simplicity of the model. 
The problem of a particle diffusing in a random poten- 
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tial, apart from being a basic problem in disordered sys- 
tems, has numerous applications in various physical pro- 
cesses including the diffusion of electrons in dis- 
ordered medium, glassy activated dynamics of disloca- 
tions in solids, dynamics of random field magnets, dy- 
namics near the helix-coil transitions in heteropolymers 
and more recently the dynamics of denaturation of a sin- 
gle DNA molecule under external force [^6|. Hence one 
expects that the knowledge of the LTD and the OTD 
in the simple Sinai type models, derived here, will have 
broader applications in several of these quenched disor- 
dered systems. 

We start with the Langevin equation of motion of an 
overdamped particle 



(2) 



where rj{t) is a thermal Gaussian white noise with zero 
mean and a correlator (ri(t)ri(t')) = 2ksT5{t — t'). For 
simplicity we set UbT = 1. The force F(x) = —dU/dx is 
derived from a potential which has a deterministic and a 
random part, U{x) = Ud{x) + U r (x). In the continuous 
version of the Sinai model we choose the random poten- 
tial U r (x) — y/a J Q ^(x')dx' to be a Brownian motion 
in space where ^(x) is a quenched Gaussian noise with 
zero mean and a correlator (^(x)^(x')} = 5(x — x'). The 
goal is to first compute the probability distribution of 
the local time, n t (a) = j Q S (x(t r ) — a) dt' and that of the 
occupation time T t {a) = f Q 9 (x(t') — a) dt' correspond- 
ing to level a for a given sample of quenched disorder 
and then obtain the disorder averaged distributions. For 
simplicity, we will consider Ud(x) to be symmetric so that 
the mean position of the particle is at zero and restrict 
ourselves to study the distributions of nt = nt(0) and 
Tt = Tt(0) corresponding to the natural choice of the 
level a = 0. However our results are easily generalizable 
to more general potentials and to arbitrary levels a. 

It turns out that the generic asymptotic scaling be- 
haviors of the LTD and the OTD, at a qualitative level, 
depend on whether the deterministic potential Ud{x) is 
unstable (Ud{x) — > — oo as x — > ±oo), stable (U^(x) — > oo 
as x — > ±oo) or flat (Ud(x) = 0). Quantitatively however, 
the LTD and the OTD do depend on the details of the po- 
tential Ud{x). To keep the discussion simple, we present 
explicit results here for the case when Ud(x) — — fx\x\, 
even though our techniques can be extended to other po- 
tentials as well. Thus we will consider Eq. (||) with the 
force F(x) — /isign(x) + y/a£_(x). For a > 0, the deter- 
ministic force is repulsive from the origin and for /i < 0, 
the force is attractive. For /i = 0, Eq. (Q) reduces exactly 
to the Sinai model. 

It is useful to start with a generalized variable 
r t = J Q V [x(t')] dt' with an arbitrary functional V [x(t')] 
which reduces to the local time n t and the occupation 
time T t when V(x) = S(x) and V(x) = 0(x) respec- 
tively. Let P x (T, t) be the probability that r t = T given 



that the particle starts at a;(0) = x. The Laplace trans- 
form Q p (x,t) = f Q °° e~ pT P x (T,t)dT plays a crucial role 
in our subsequent analysis. By definition Q p {x,t) — 

(e p Xi v [ x ^ )] dt \ where () denotes the average over all 
histories of the particle upto time t starting at x at t = 0. 
Using the evolution equation (g) it is straightforward to 
see that Q p (x,t) satisfies, for arbitrary F(x), the back- 
ward Fokker-Planck equation 



OQp = l d 2 Q p 
dt 2 dx 2 



F(x) 



dQ p 
dx 



pV(x)Q p 



(3) 



with the initial condition Q p (x,0) = 1. After a fur- 
ther Laplace transform, now with respect to t, u(x) — 
Io° G~ at Q p {x,t)dt satisfies the equation 

1 



-it + F(x) u — [a + pV(x)] u 



-1. 



(4) 



where u (x) = du/dx and we have suppressed the a and 
p dependence of u(x) for notational convenience. So far 
the discussion is quite general. We now consider the local 
and the occupation time separately. 

Local Time : In this case V(x) — 5(x). We need 
to solve Eq. (||) separately for x > and x < 
and then match the solutions at i = 0. We write 
u±(x) — l/a + A±y±(x) where y±(x) satisfy the ho- 
mogeneous equations 



1 „ 

- 2 v± 



F (x)y± - a v± = o, 



(•5) 



respectively in the regions x > and x < with the 
boundary conditions y+{x — * oo) = and y~{x — > 
— oo) = 0. The constants A± are determined from the 
matching conditions, u + (0) — u_(0) = u(0) and u + (0) — 
u_(0) — 2pu(Q). Eliminating the constants, we get 
w(0) = X(a)/a\p + X(a)] where A (a) = [z_(0) -z+(0)]/2 
and z±(x) — y ± {x)/y±{x). For simplicity, we will re- 
strict ourselves only to Po(n t = T,t) corresponding to 
the natural choice of the starting point x(0) — 0, though 
our methods can be easily generalized to arbitrary initial 
positions. Since A(a) is independent of p, one can easily 
invert the Laplace transform u(0) with respect to p to 
get 



G(a) - / , 

10 



' P (rit =T,t)dt 



A(a) 



-A(a)T 



. (6) 



a general result valid for arbitrary F(x) . To proceed more 
we choose F{x) — /isign(x) + y/a^x) and consider the 
implications of Eq. (||) for the pure case a = first for 
arbitrary a. 

Solving Eqs. (Q) for F(x) = /isign(x) we find y±(x) = 
y±(Q)e Tl3x with (3 = n + vV + 2a. This gives A(a) = 
(3 = fi+ \J + 2a. By analysing this Laplace transform 
one finds that there are three different asymptotic (large 
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t) behaviors of Pq(T, t) depending on whether the poten- 
tial is unstable (/i > 0), stable (p, < 0) or flat (// = 0). 
These asymptotic behaviors for the particular potential 
chosen here (Ud(x) = p\x\) can be shown to be typical 
for generic potentials. 

• Unstable potential (/i > 0) : In this case, the LTD 
approaches a steady state in the large window size t — > oo 
limit, P (T) = 2^e~ 2 ^ T obtained by taking a 
limit in Eq. (0). Physically it indicates that for re- 
pulsive force the particle eventually goes to either oo or 
—oo and occasionally hits the origin according to a Pois- 
son process. This asymptotic exponential distribution, 
Pq(T) = A(0)e -A (°) T is indeed universal (up to a rescal- 
ing factor of time) for any unstable potential. 

• Stable potential (// < 0) : For generic stable po- 
tentials the system approaches a stationary state in the 
large t limit and the stationary probability distribution 
p(x) for the position of the particle is given by the Gibbs 
measure, p(x) — e~ Wi ^ x > /Z where Z = e^ 2Ud ^dx 
is the partition function. Hence as t — > oo, simple er- 
godicity arguments indicate that the local time n t — > 
f*(5[x(t')])dt' -+p(0)t, i.e., the LTD is simply P (n t = 
T,t) = S (T — p(0)t), a result that can also be proved rig- 
orously fL7| . For example, for the example chosen here, 
one finds p(0) = \p\ and hence Po(T,t) = S(T — \p\t) as 
t — > oo. It turns out that the approach to this asymp- 
totic delta function is rather slow in the following sense. 
Indeed one finds that in the scaling limit t — > oo, n t — > oo 
but keeping the ratio n t /t — y fixed, the LTD behaves as, 
Po{ n t =T,t) ~ exp [— i$(T/t)] where the large deviation 
function (LDF) $(y) is given from Eq. ^ via a steepest 
descent calculation 



$(y) = max Q [-a + y\{a)} , 



(7) 



which, once again, is a general result valid for any ar- 
bitrary stable potential. In our example, substituting 
A (a) = — \p\ + y/[i 2 + 2a in Eq. (Q) and maximizing we 
get the exact LDF, = (y- \fi\) 2 /2. Thus in this par- 
ticular example Po(T,t) is a shifted Gaussian with mean 
at T/t = \p\ whose width decreases as a power law i -1 / 2 
and only in the strict t — * oo limit, it approaches a delta 
function peaked at T/t = \p\. This result also follows 
if one assumes the validity of the central limit theorem 
(CLT). However, for arbitrary stable potential, the CLT 
breaks down and the LDF $(1/) shows departure from 
the simple quadratic form |l7| , |l8| ]. 

• Flat potential (/i = 0) : In this case, one finds by 
inverting the Laplace tranform in Eq. (||) that the LTD 
is Gaussian for all T and t, P (T,t) = y/2/ni e~ T * '/ 2t . 

These behaviors for the pure system [a — 0) get dras- 
tically modified when the random potential is switched 
on (<7 > 0). The equation (^J) still remains valid for 
each realization of F(x). For instance by taking the 
a — > limit one may recover the expression of the ther- 
mal average of the local time which was interpreted in 



]15| as a trapping time. Here our aim is quite different 
since we want to compute the disorder averaged LTD 
Po(n t =T,t). From Eq. (0), one needs to know the dis- 
tribution of X(a) — [z-(0) — z_|_(0)]/2 which is now a 
random variable since F (x) is random. It turns out the 
distribution of A (a) can be computed exactly by adopting 
techniques that have appeared before in the Sinai model 
in other contexts [ p~5|J]~S| , |20| ] . We defer the technical de- 
tails for a future publicati on ]l7[ and p resent only the 
final results here. One gets exp[— X(a)T] — q 2 (T) with 



K 



q(T) = (l+aT)-^ 2a - 



n/a 



K 



p. I a 



(8) 



where K v (x) is the modified Bessel function of order v 
ifjlf . Averaging Eq. (||) over disorder we finally get the 
exact formula 



P (T,t)e- at dt 



1 d r o , . , 
a dl 



(9) 



where q(T) is given by Eq. (||). The asymptotic behav- 
iors can be deduced by analysing Eq. (|^). 

• fx > : In this case, we find that as t — > oo, 
P (T,t) tends to a steady state distribution, Pq{T) — 
2(1(1 + aT)- 2 ^/"- 1 for all T > 0. Thus the disorder 
averaged LTD has a broad power law distribution even 
though for each sample the LTD has a narrow exponen- 
tial distribution. This indicates wide sample to sample 
fluctuations and lack of self averaging. 

• (j, < : In this case, for each sample of the dis- 
order, the LTD tends to the delta function Pq{T, t) — > 
6 (T — p(0)t) as discussed before. However, the Gibbs 
measure p(0) varies from sample to sample. On aver- 
aging over disorder (or equivalently the peak positions), 
one finds a broad distribution for Pq (T, t) . In the scaling 
limit t — ► oo, n t = T — * oo but keeping the ratio T/t 
fixed, we find that Po(T,t) — > jf(T/t) where the scaling 
function can be computed exactly by analysing Eq. (^), 



f(y) 



( j(2 ( 7 3 )(^ 1 )/ 2 r 2 (i/) 



y^-^e-y^W^y/a), 
(10) 



with v = \n\/a and W v ,u(x) is the Whittaker function 
pl[ . The scaling function increases as f(y) ~ y v ~ x for 
small y and eventually decays for large y as f(y) ~ 
y(4i/-3)/2 e -2t//er_ Once again the disorder modifies the 
behavior of the LTD rather drastically. 

• /i = (Sinai model) : In this case we find that for 
large t, P (n t =T,t) -> ^- t f s (T/t) where f s (y) = 

2e~ y / a Ko(y /a)/y. However this scaling breaks down for 
very small y when y << a. 

We now turn to the OTD, R (T,t) = Prob(T t = 
T,t) given that the particle starts at i = at 
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t = 0. In this case, the double Laplace transform 
Up (x) = J °° dte~ at J*e-P T R x (T,t)dT satisfies Eq.§) 
with I^(cc) = 8(x). Since the rest of the calculations are 
very similar to the LTD case, we just present the final 
results omitting the details. 

• Unstable potential: Consider the pure case (a = 0) 
first. Since the deterministic potential is symmetric, one 
has R Q (T,t) = Ro(t- T,i), i.e. the OTD is symmet- 
ric around its mean value t/2. In the limit of a large 
window size t — > oo, it turns out that the part of the 
OTD to the left of the midpoint T = t/2 approaches 
a steady (t independent) distribution Rr,(T) with the 
normalization f™ R L {T)dT = 1/2. The right half of 
the OTD, which carries an equal total weight 1/2 is 
pushed to oo since the midpoint t/2 itself goes to oo. 
This conclusion is valid for any symmetric deterministic 
potential. For the case F{x) = /isign(x), we get ex- 



plicitly, R L (T) 



9 _ 

// e 



1 — 3y/nue 9u erfc(3u) / ^/ttu 



where u = //yT/2 and erfc(x) is the complementary er- 
ror function. Thus R L (T) w n^2/irT for small T and de- 
cays exponentially for large T, R L (T) ~ t^I 2 e~^ T I 2 . 
When the disorder is switched on [a > 0), this asymp- 
totic behavior for the pure case does not change qualita- 
tively. Once again the left part of the disordered averaged 
OTD tends to a t independent form Rt,(T). In fact, the 
small T behavior of RjT) ps ^2/ttT remains the same 
as in the pure case. However, for large T, while the OTD 
still decays exponentially Rl{T) ~ e~ bT , the decay coef- 
ficient b turns out to be different from the pure case and 
is given by the negative of the real part of the zero (clos- 
est to the origin) of K fl / cr {^/2z/a) = in the complex z 
plane. 

• Stable potential: As in the case of the LTD we find 
that for the pure case, for generic stable potential Ud(x), 
the OTD approaches a delta function in the t — > oo limit, 

R Q (T t = T,t) 6 (T - ^-tj where Z is the equilib- 
rium partition function and Z + = e^ 2Ud ^dx is the 
restricted partition function. This result again follows 
from simple ergodicity arguments. Also, the approach 
to this asymptotic form is slow and one finds the scaling 
behavior, R {T t = T,t) ~ e ~**( T /*) in the scaling limit 
T — > oo, t — ► oo keeping T/t fixed where ^(y) is the 
LDF. For the case F(x) = /isign(a;), the LDF can be 
computed exactly and we get #(y) = 2fi 2 (y — 1/2) 2 . Re- 
cently the LDF has been computed exactly for other sta- 
ble potentials also pS| . In presence of disorder (cr > 0), 
this asymptotic behavior gets modified as in the case of 
LTD and we find R (T,t) w \ f (T/t) in the scaling limit. 
The exact calculation of the scaling function f (T/t) is 
nontrivia! but the final answer turns out to be a decep- 
tively simple Beta law, 



where v = \fi\/cr and B(v : v) is the standard Beta func- 
tion . If one tunes the parameter v by either varying 
H or the disorder strength a, this OTD exhibits an inter- 
esting phase transition in the ergodicity of the particle 
position at v c = 1. For v < v c , the distribution in Eq. 
( pi] ) is concave with a minimum at y = 1/2 and diverges 
at the two ends y = 0,1. This means that paths with 
small number of zero crossings (such that T is close to 
either or t) carry more weights than the paths that cross 
many times (for which T is close to t/2), i.e. the particle 
tends to stay on one side of the origin as in the case of a 
Brownian motion. Exactly the opposite situation occurs 
for v > v c where f (y) is maximum at its mean value 
1/2 indicating largest weights for paths that spend equal 
times on both sides of x = 0. It is interesting to notice 
that similar types of Beta laws also arise in the study of 
certain perturbed Brownian motion ]2^] . 

• Flat potential (fi = 0) : For the pure case (<r = 0), our 
method reproduces the 'arcsine' law for the OTD of an 
ordinary Brownian motion, Ro(T,t) = l/TTy/T(t — T). 
In the presence of disorder (cr > 0), i.e. for the the Sinai 
model, we find that the left part of the OTD for < T < 
t/2 has the large t behavior, R%(T,t) w ^R(T). The 
right half of the OTD for t/2 < T < t is just the symmet- 
ric reflection of the left part. The t independent function 
R(T) has a complicated form but with simple limiting 
behaviors, R{T) » a^2/nT as T -> and R(T) w 1/2T 
for large T, consistent with the normalization condition 
Jl' 2 R%{T,i)dT = 1/2. 

In summary we have studied analytically the LTD and 
OTD for a particle moving in an arbitrary one dimen- 
sional potential and shown how the quenched disorder 
changes the asymptotic behaviors of these distributions. 
Our exact results are consistent with the general notion 
that 'disorder broadens distributions' of physical quanti- 
ties. Recently several asymptotic exact results for other 
quantities in the Sinai type models were derived using 
a real space renormalization group (RG) treatment [p3| . 
Reproducing the exact results presented here either via 
the RG method or by the replica method and extending 
our results to higher dimensions remain as challenging 
open problems. 

We thank M. Yor for useful discussions. 



foiv) 



B {u, v) 



W-y)] 1 



0<y<l (11) 
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